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Abstract. 

With the aid of simple numerical models, we discuss a particular aspect of the 
interaction between stellar orbital periods inside elliptical galaxies (Es) and the 
parent cluster tidal field (CTF), i.e., the possibility that collisionless stellar evapo- 
ration from Es is an effective mechanism for the production of the recently discov- 
ered intracluster stellar populations (ISP). These very preliminary investigations, 
based on idealized galaxy density profiles (such as Ferrers density distributions) 
show that, over an Hubble time, the amount of stars lost by a representative 
galaxy may sum up to the 10% of the initial galaxy mass, a fraction in interesting 
agreement with observational data. The effectiveness of this mechanism is due 
to the fact that the galaxy oscillation periods near its equilibrium configurations 
in the CTF are of the same order of stellar orbital times in the external galaxy 
regions. 
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1. Introduction 

Observational evidences of an Intracluster 
Stellar Population (ISP) are mainly based 
on the identification of intergalactic plan- 
etary nebulae and red giant branch stars 
(see, e.g., Theuns & Warren 1996, Mendez 
et al. 1998, Feldmeier et al. 1998, Arnaboldi 
et al. 2002, Durrell et al. 2002). Overall, 
the data suggest that approximately 10% 
(or even more) of the stellar mass of the 
cluster is contributed by the ISP (see, e.g., 
Ferguson & Tanvir 1998). 



The usual scenario assumed to ex- 
plain the finding above is that gravita- 
tional interactions between galaxies in clus- 
ter and/or interactions between the galax- 
ies with the tidal gravitational field of the 
parent cluster (CTF), lead to a substantial 
stripping of stars from the galaxies them- 
selves. 

Here, supported by a curious coinci- 
dence, namely by the fact that the char- 
acteristic times of oscillation of a galaxy 
around its equilibrium position in the CTF 
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are of the same order of magnitude of the 
stellar orbital periods in the external part of 
the galaxy itself we suggest tliat an addi- 
tional effect is at work, i.e. we discuss about 
the possible "resonant" interaction between 
stellar orbits inside the galaxies and the 
CTF. 

In fact, based on the observational 
evidence that the major axis of cluster 
Es seems to be preferentially oriented to- 
ward the cluster center, N-body simula- 
tions showed that galaxies tend to align re- 
acting to the CTF as rigid (Ciotti & Dutta 
1994). By assuming this idealized scenario, 
a stability analysis then shows that this 
configuration is of equilibrium, and allows 
to calculate the oscillation periods in the 
linearized regime (Ciotti & Giampieri 1998, 
hereafter CG98). 

Prompted by these observational and 
theoretical considerations, in our numerical 
explorations we assume that the galaxy is 
a triaxial ellipsoid with its center of mass 
at rest at the center of a triaxial cluster. 
The considerably more complicate case of 
a galaxy with the center of mass in rotation 
around the center of a spherical cluster will 
be discussed elsewhere (Ciotti & Muccione 
2003, hereafter CM03). 



2. The physical background 

Without loss of generality we assume that 
in the (inertial) Cartesian coordinate sys- 
tem C centered on the cluster center, the 
CTF tensor T is in diagonal form, with 
components Ti (i = 1, 2, 3). By using three 
successive, counterclockwise rotations {tp 
around x axis, i? around y' axis and iJj 
around z" axis) , the linearized equations of 
motion for the galaxy near the equilibrium 
configurations can be written as 



.. AT32A/32 
^ = J V, 

1!) = —1), 



h 

AT21A/2 



where AT is the antisymmetric tensor of 
components ATy = Ti — Tj , and li are the 
principal components of the galaxy inertia 
tensor. In addition, let us also assume that 
Ti > T2 > T3 and h < h < h, i.e., that 
AT32,Ar3i and Ar2i are all less or equal 
to zero (see Section 3). Thus, the equilib- 
rium position of is linearly stable, and its 
solution is 



= ipM COs(tJ^t), 
d = COs{uJi)t), 
Tp ^ TpM COS^LO^t), 

where 
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For computational reasons the best ref- 
erence system in which calculate stellar or- 
bits is the (non inertial) reference system 
C" in which the galaxy is at rest, and its 
inertia tensor is in diagonal form. The equa- 
tion of the motion for a star in C is 

a;' = n^x-2nAv'-nAx'-nA{nAx'), (4) 

where x — Tl{(p,'d,ip)x' , and 

n = ((yicostJcos + i^sin?/;, — i^cosiJsinV' 
+d cos ^p,(p sin 1} + -tp). (5) 

In eq. (4) 



n'^x 



(7^^T7^)a;^ 



(6) 



(1) 



where <pg{x') is the galactic gravitational 
potential, Vx' is the gradient operator in 
C", and we used the tidal approximation 
to obtain the star acceleration due to the 
cluster gravitational field. 

3. Period estimations 

For simplicity, in the following estimates we 
assume that the galaxy and cluster densi- 
ties are stratified on homeoids. In partic- 
ular, we use a galaxy density profile be- 
longings to the ellipsoidal generalization of 
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the widely used 7-models (Dehnen 1993, 
Tremaine et al. 1994): 



Pg(m) 



3-7 



(7) 



where Mg is the total mass of the galaxy, 
< 7 < 3, and 



m 



E 

i=l 



at 



ct\>a2> 0:3. (8) 



The inertia tensor components for this 
family are given by 



All 
"3 



(9) 



where hg = Pg{m)m'^dm, and so /i < 
I2 < h- Note that, from eq. (3), the fre- 
quencies for homeoidal stratifications do 
not depend on the specific density distri- 
bution assumed, but only on the quantities 
(ai,a2,Q^3)- We also introduce the two el- 
lipticities 



^ = 1 
ai ~ 



^ = 1 

ai ~ 



V, 



(10) 



where c < r] < 0.7 in order to reproduce 
realistic flattenings. 

A rough estimate of characteristic stel- 
lar orbital times inside m is given by 
Porb{m) ~ 4Pdyn(m) = 3tt / Gp^{m) , 
where 'p^{m) is the mean galaxy density in- 
side m. We thus obtain 



„ /a?i(l-e)(l-r;) 

xmT/2(l + m)(3-'')/2 yrs, (11) 

where Mg^n is the galaxy mass normalized 
to IO^^Mq, ai,i is the galaxy "core" ma- 
jor axis in kpc units (for the spherically 
symmetric 7 = 1 Hernquist 1990 models, 
i?e — l.Sai); thus, in the outskirts of nor- 
mal galaxies orbital times well exceeds 10^ 
or even 10^ yrs. 

For the cluster density profile we as- 
sume 



/ \ Pc.a 

(1 + m'^Y 



(12) 



where m is given by an identity similar to 
eq. (8), with ai > 02 > 03, and, in analogy 
with eq. (10), we define 02/01 = 1 — /i and 
os/ai = l — u, with fi< f <1. 

It can be shown that the CTF com- 
ponents at the center of a non-singular 
homeoidal distribution are given by 



-2TTGpc,0Wi{p-, v) 



(13) 



where the dimensionless quantities Wi are 
independent of the specific density profile, 
Wi < «'2 < for oi > 02 > 03, and so the 
conditions for stable equilibrium in eq. (1) 
are fulfilled (CG98, CM03). 

The quantity pcfl is not a well mea- 
sured quantity in real clusters, and for its 
determination we use the virial theorem, 
Mc<Jy = —U, where cry is the virial ve- 
locity dispersion, that we assume to be es- 
timated by the observed velocity disper- 
sion of galaxies in the cluster. Thus, we can 
now compare the galactic oscillation peri- 
ods (for small galaxy and cluster flatten- 
ings, see CM03) 



= 

P« = 
P^ = 
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(where 01,250 = ai/250kpc, and try, 1000 = 
(Tv/lO^km s^^) with the characteristic or- 
bital times in galaxies. From eqs. (11) and 
(14), it follows that in the outer halo of 
giant Es the stellar orbital times can be 
of the same order of magnitude as the os- 
cillatory periods of the galaxies themselves 
in the CTF. For example, in a relatively 
small galaxy of Mg,ii = 0.1 and ai,i = 1, 
Porb — 1 Gyr is at m ~ 10 (i.e., at ~ 5i?c), 
while for a galaxy with Mg,ii = 1 and 
ai,i = 3 the same orbital time is at m ~ 7 
(i.e., at ~ 3.5Re). 

4. The numerical integration scheme 

Here we present the preliminary results 
obtained by using only a very idealized 
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Fig. 1. Histogram of m at t vs. A'osc for a galaxy of 7\fg = 10^" Mq and ai = 20 
kpc, with e = 0.2 and rj — 0.4. Tlie cluster parameters are = 0.2, v = 0.4, ai,250 — 
cv,iooo — 1, and the galaxy maximum oscillations angles in eqs. (2) are fixed to 0.2 rad. 
In this simulation iVtot — 10^. 



galactic density profile, while more realis- 
tic galaxy density distributions, such as tri- 
axial Hernquist models, are described else- 
where (CMOS, Muccione & Ciotti 2003, 
hereafter MC03). In particular, here we ex- 
plore the evaporation process in the case 
of a Ferrers (see, e.g., Binney & Tremaine 
1987) models, with density profiles given by 



Pa 



p<,(0)(l-m2)« form<l 

for TO > 1 ^ ' 



where to is the homeoidal radius defined 
as in eq. (8). The case n — corresponds 
to an anisotropic harmonic oscillator, while 
for n > the density distribution is radi- 
ally decreasing. A nice property of Ferrers 
models (for integer n) is that their gravi- 
tational potential can be simply expressed 
in algebraic form (see, e.g., Chandrasekhar 
1969). 

To integrate the second order differ- 
ential equations (11) for each star, we 
use a code based on an adapted Runge- 
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Kutta routine. The initial conditions are 
obtained by using the Von Neumann re- 
jection method, and with this method we 
arrange iVtot initial conditions in configu- 
ration space which reproduce the density 
profile in eq. (15). In MC03 and CMOS the 
three components of the initial velocity for 
each star are assigned by considering the lo- 
cal velocity dispersion of the galaxy model 
at rest, while here, for simplicity, each star 
at t = is characterized by null velocity. 
Within t = in (whore t^ = 1.5 x 10^" 
yrs), the code checks which of the initial 
conditions result in a orbit with m > 1. 
This escape condition is very crude, and 
a more sophisticated criterium is adopted 
in MC03 and CMOS, where untruncated 
galaxy models arc used. In standard sim- 
ulations we use, as a rule, iVtot = 10^, thus, 
from this point of view, we are exploring 
10"'' independent "1-body problems" in a 
time-dependent force field. 

The simulations are per- 
formed on the GRAVITOR, the 
Geneva Observatory 132 processors 
Beowulf cluster (http://obswww. 
unige . ch/ ~ pfennigs/ gravitor / 
gravitor_e.html) . 

5. Preliminary results and conclusions 

In Fig. 1 we show the results of one of our 
preliminary simulations for a galaxy model 
with n = 1, Mg = 10^^ Mq, semi-mayor 
axis ai = 20 kpc, flattenings e = 0.2, 
T] = 0.4, and maximum oscillation angles 
equals to 0.2 rad. The cluster parameters 
are ai,250 = o-vaooo 1, /x = 0.2, v = 0.4, 
and the total number of explored orbits is 
^tot = 10^. In the ordinate axis we plot the 
number of escapers as a function of their 
homeoidal radius at t = 0. Note how the 
zone of maximum escape is near m ~ 0.8, 
and how the total number of escapers is a 
significant fraction of the total number of 
explored orbits (actually, with the galaxy 
and cluster parameters adopted in the sim- 
ulations, A^osc — O.lA^tot)- This number is a 
somewhat upper limit of the expected num- 
ber in more realistic simulations: in fact, 



1) as described in Section 4, we adopted 
a very weak escape criterium, and 2), in 
more realistic galaxy models the density de- 
crease in the outer parts is stronger than 
in Ferrers models, and so we expect that a 
smaller number of stars will be significantly 
affected by the CTF. Both points are ad- 
dressed and discussed elsewhere, in MCOS 
in a preliminary and qualitative way for an 
Hernquist model at the center of a clus- 
ter, while in CMOS we present the results 
of a systematic exploration of the parame- 
ter space for untruncated galaxies both at 
the center and with the center of mass in 
circular orbit in a spherical cluster. 

In any case, the simple model here ex- 
plored looks promising, with a number of 
escapers in nice agreement with the obser- 
vational estimates. 
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